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P. ■ Abstract 



PrBa2Cu307 (PBCO) is unique in the RBa2Cu307 (R=rare earth) series 
because it is not superconducting. In fact, for Yi_a;Pra;Ba2Cu307, Tc drops 
monotonically with Pr concentration, with Tc going to zero at ~55% Pr. 
There have been many studies of this material with the hope that an expla- 
nation for the lack of superconductivity in PBCO might help explain why 
YBa2Cu307 (YBCO) is superconducting. Explanations center around hole 
localization, requiring an extra hole on the Y (Pr) site or a localized hole in 
the O 2p shell, at the expense of a mobile hole in the Cu-0 planes. To help 
provide clues that could point to a particular model, and to search for anoma- 
lies in the local structure, we present K-edge XAFS (X-Ray Absorption Fine 
Structure) data for various concentrations of Pr in Yi_a;Pr^Ba2Cu307. The 
character of the Pr K-edge XAFS data indicate that most of the Pr substitutes 
onto the Y site and is well ordered with respect to the unit cell. These data 
also show that the amplitude of the first Pr-0 peak is greatly reduced when 
compared to the first Y-0 peak in pure YBCO, and decreases with increasing 



Pr concentration. In contrast, the Y K-edge data for these alloys show little 
if any change in the oxygen environment, while the Cu K-edge data show a 
10% reduction in the first Cu-0 peak. Fits to the Pr data suggest that some 
oxygen atoms about the Pr become disordered and/or distorted; most of the 
Pr-0 nearest-neighbor distances are 2.45 A, but about 15-40% of them are in 
a possibly broadened peak at 2.27^0^2 A. The Cu K-edge XAFS data show 
a slight broadening but no loss of oxygens, which is consistent with a radial 
distortion of the Pr-0 bond. The existence and the size of these two bond 
lengths is consistent with a mixture of Pr^+ and Pr^"*" bonds, and to a formal 
valence of +3.33^Q']^g for the Pr ion. 

PACS numbers: 74.70.Vy, 74.75.+t, 78.70.Dm, 61.70.-r 
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I. INTRODUCTION 

An interesting property of YBa2Cu307 (YBCO) is that one may replace yttrium with 
almost any rare earth and still obtain a high-T^ superconductor Jllo isostructural with YBCO. 
However, there are two rare earths that do not form a high-T^ superconductor when attempts 
are made to substitute them onto the Y site: cerium and praseodymium.l§Q Attempts to 
substitute with Ce have not produced the YBCO structure, even in small concentrations. 
Praseodymium, however, substitutes readily onto the Y site, with Tc dropping monotonically 
with Pr concentration. Therefore, PrBa2Cu307 (PBCO) is a unique material in that it has 
the YBCO structure, with a rare earth instead of Y, yet is not superconducting. In contrast 
substitution by Gd, for example, has essentially no change in Tc or in the crystal structure 
There has been a large effort to explain this reduction in Tc as a new avenue in the search 
to explain why YBCO is superconducting in the first place. PBCO is interesting in its own 
right, however, as an insulating layer between layers of YBCO in multilayer films, or as a 
well matched substrate for epitaxial films of YBCO. 

There are many related explanations for the decrease in Tc in Yi_^Pr^.Ba2Cu307 
(YBCO:Pr), mostly involving the localization of holes through valence change or magnetic 
pair breaking mechanisms. The earliest arguments note that Ce, Pr and Tb are the only 
rare earths that can exist in either a +3 or +4 valent state.a Since Y only exists in a +3 
valent state, one explanation for the decrease in Tc with Pr concentration is that some 
fraction of the Pr exists in its formal +4 state, filling a hole state in the Cu-0 planes and 
localizing it onto the Pr site. There is some experimental evidence to support this view, 
notably the measured magnetic susceptibility of 2.7 fiB per Pr ion in PBCO,Q which may be 
compared to the Pr^+ ion (3.58 /i^) and the Pr"^"*" ion (2.54 /x^). However, this result should 
be considered in light of some model systems that give moments significantly different than 
the ionic moments, such as the moment of BaPrOa (0.7 /i_B),l3 where Pr"^"*" is presumably the 
natural state. 

Another possible mechanism for suppressing Tc in this material utilizes the fact that, 
except for Ce, the Pr^+ ion is the largest of all the rare earth +3 ions. The 4/ electrons are far 
enough extended to hybridize with the O 2p electrons, and thus interfere with the conduction 
band.u In fact, crystal field measurements can be reconciled with Pr in a completely +3 
state,&u with broad transitions that can be attributed to hybridization between the Cu 5d-0 
2p states and the Pr 4/ states.Ll 

X-Ray Absorption Near-Edge Structure studies (XANES) also indicate that Pr is pri- 
marily trivalent. Comparisons of the Pr L///-edge of PBCO to that for Pr02 and Pr203 
indicate a valence close to +3,t3 or slightly above it.lllHl3 The latter studies suggest that a 
small amount of Pr +4 is present; Horn et alM give a mean valence of -1-3.1, while Lytle et 
alM find that the valence decreases from +3.45 for YBCO:Pr 20% to +3.25 for YBCO:Pr 
60% with all of the Pr^+ on the Cu(2) site (see Fig. |I] for structural notation). M-edge 
studiesB and electron-energy loss spectroscopy (EELS) studiealil also find that the Pr va- 
lence is close to +3. These same studies show that the number of holes in the CuO? planes 
does not change with Pr concentration, while normal-state Hall effect measurement al3 show 
a decrease in the number of mobile charge carriers with x, indicating that the holes are be- 
ing localized rather than filled. Photoemission experiments also find that the Pr valence is 
close to +3 and that there is significant hybridization of the Pr 4/ states with other valence 



states.El In addition, thin films of Yo.5Tbo.5Ba2Cu307 liave recently been produced, with a 
Tc of 92 K.EJ Such a high T^ for a 50% terbium sample brings into question the argument 
that the suppression of T^ arises solely because two valence states can exist. 

A theory which combines many of these ideas, proposed by Fehrenbacher and Rice,ll^ 
shows that three electronic configurations are favored: a Pr'^^(4/)-0 state, a Pr^"'"(4/)-0 
state, and a Pr^"''(4/L)-0, where the Pr'^+(4/L) indicates that the hole does not reside 
within the praseodymium's electronic structure; rather, the hole is in the ligand. The low- 
energy ligand hole is obtained via a superposition of O 2pt^ orbitals with / symmetry. The 
most energetically favorable configuration has Pr^+ 60-70% of the time, with the remaining 
Pr-0 bonds split between the other two configurations. They reason that such a mixture 
would be consistent with all the above measurements, since ionic Pr^+ would only exist 
15-20% of the time, even though formally the Pr is in a -|-4 state 30-40% of the time. 

There is also some structural evidence from neutron and x-ray diffraction for a mixed Pr 
valency. Diffraction results show that PBCO is isostructural with YBCO. The diffraction 
result for the average Pr-0 nearest-neighbor distance of 2.45 AcHc3 is reasonably consistent 
with Pr'^+-0 bond lengths in other Pr oxides, such as in bixbyite Pr203, which has a range of 
Pr-0 bond lengths from 2.33 A to 2.66 A.IH3 It has been noted, however, that both the Cu-0 
plane separationc3'E3 and the Pr-0 (2) and Cu(2)-0(3) bond lengthsc3 are more consistent 
with a valence in the range of +3.3-3.4, as determined by interpolating trends in these bond 
lengths with rare-earth ionic radius. However, this interpolation argument has at least two 
possible flaws. First, just because the trend with ionic radius is no longer linear for Pr does 
not necessarily mean Pr is in a mixed valence state. The Pr^"*" ion is the largest ion anyone 
has been able to place between the Cu-0 planes in this structure, so new forces, including 
more oxygens bonding to the Pr than, say, Y, may be coming into play. In addition, there 
are several compounds that show a wide range of Pr^+ bond lengths, such as bixbyite Pr203. 
Secondly, if three Pr-0 states exist, as suggested by Fehrenbacher and Rice,ll3 then one would 
expect to see two, and perhaps three, different Pr-0 bond lengths; Pr^+-0 bonds range in 
length from 2.18 A in BaPrOa to 2.32 A in Pr02. Therefore the three Pr-0 states should 
have bond lengths made up of a long Pr^+(4/)-0, a short Pr^+(4/)-0, and a Pr^+(4/L)-0 
(formally, Pr^"'"-0), whose bond length should be closer to Pr^"*"(4/)-0. Diffraction would 
not be sensitive to such a distribution. It will either see the dominant bond length (perhaps 
the Pr^+(4/)-0), or an average bond length consistent with a mixed valency, depending on 
the degree of disorder. 

Likewise, trends in the local O environment around the Pr with Pr concentration are 
only obtainable with large Pr concentrations. If trends exist in the concentration region 
where YBCO:Pr is still superconducting, that is, Pr concentrations < 55%, diffraction may 
not see them. Therefore diffraction results do not rule out disorder and/or distortions of the 
O around the Pr. A distortion in the Pr local structure could indicate the presence of some 
formal or ionic Pr^+, but other disorder may also exist. For instance, there is experimental 
evidence from x-ray diffraction to show that the Cu(2)-0(4) distance is anomalously short,lHll 
and therefore may impede charge transfer between the planes and the chains. The length of 
this bond is consistent with the trend with ionic radius of Pr^+jB^ but it is still the shortest 
Cu(2)-0(4) bond measured in the YBCO-type materials. Any of these results could produce 
some localization of charge carriers. 

Since local distortions and/or disorder in the oxygen environment around praseodymium 



may indicate the presence of some Pr with a formal valence of 4+, we have con- 
ducted X-Ray Absorption Fine Structure (XAFS) experiments on several concentrations 
of Yi.^Pr^BaaCuaOr, at the Cu, Y, and Pr K-edges (E^ ~ 9 keV, 17 keV, and 42 keV, 
respectively). (A preliminary version of this work has recently appeared in the literature.Bj) 
XAFS is ideally suited to this problem, because it can give us precise information about the 
local structure around the absorbing atom, including average bond distances and harmonic 
broadening factors. In particular, we should be able to determine exactly which site the Pr 
is occupying and observe changes in the XAFS compared to YBCO. By taking data at three 
different absorption edges, we can further constrain our results for the Cu02 planes and get 
more information regarding the 0(4) site. 

In Sec. H we describe some details of the sample preparation and data collection. We 
explain our data reduction procedures and describe some visible features of the data in Sec. 
p| . Our analysis and fit results are detailed in Sec. |^ We discuss the implications of our 



results in Sec. |V] and give a brief conclusion of our findings in Sec. |VI|. 

II. SAMPLES AND DATA COLLECTION 
A. XAFS Samples 



All the samples were prepared by the same procedure as described in Ref. ^. Stoi- 
chiometric amounts of CuO, BaCOs, Y2O3, and PreOn were mixed, ground, and heated 
at 920°C in air for 20 hours. They were then reground, pressed into pellets, calcined at 
940°C in an oxygen atmosphere, and slowly cooled at 60°C/h to room temperature. All the 
Yi_a,Pra,Ba2Cu307_y samples were characterized by x-ray diffraction using Cu Ka radiation. 
The patterns indicate that all the samples with x < 0.6 are essentially single phase. The 
100% Pr sample may have small amounts of phase impurities, but they are too small to be 
seen in the XAFS (see Sec. P^ ). Oxygen content was analyzed by an improved gaseous vol- 
umetric technique and found to be within the range of 6.95±0.02 for all samples. In all the 
compounds T^ was determined by a standard four-probe method. Superconducting phase 
purity was checked using an AC mutual inductance bridge. Only one step was observed for 
all superconducting samples. Fig. ^ shows Tc as a function of Pr concentration and Table | 
gives the lattice parameters and the oxygen content measurements for these samples. 

XAFS scans were performed on the 0, 30, 50, and 100% Pr samples at about 80 K for 
the Pr K-edge and about 50 K for the Cu and Y K-edges. To prepare the samples for these 
transmission experiments, we reground the pellets and sifted the powder through a 30 /xm 
filter onto pieces of scotch tape. The tape was then cut and stacked into multiple layers 
to obtain samples of thickness, t, such that /it ~ 1 at the absorption edge, where fi is the 
absorption coefficient. All XAFS experiments were performed at SSRL on beamlines 7-3 
and 10-2. Measurements of the incident and transmitted intensities were made with gas- 
ionization chambers. Si (111) monochromator crystals were used for the Cu and Y K-edge 
experiments, while Si (400) crystals were required for the Pr K-edge. 

III. GENERAL FEATURES OF THE DATA 



A. Standard reduction and comparison of different absorption edges 

Y, Pr, and Cu K-edge data for x(k) = (/i - l^o)/ l^o show XAFS oscillations beyond 16 
A~^ for all K-edges. Fig. ^ shows kx(k) for all three K-edges for Yo,5Pro.5Ba2Cu307. The 
"free atom absorption" /Xq is determined in the usual wayn3 by fitting a spline function or a 
polynomial through the data above the absorption edge. The photoelectron wave vector k is 

given by k= \J {2'm^/h'^){E — Eo), where Eq is the K-edge threshold energy. This procedure 
was modified somewhat for the Pr K-edge data, as detailed in part B of this section. 

The most striking difference between the Y and the Pr K-edge XAFS in Fig. |^ is the 
overall drop in amplitude of the Pr K-edge XAFS, particularly in the low k part of the 
spectrum. The amplitudes of the high-k XAFS are similar. The Cu K-edge XAFS are 
comparable to data we have previously obtained on YBCO. 

The Fourier transform (FT) of the data yields peaks in r-space corresponding to different 
radial distances from the excited atom. The peaks in the FT of kx(k) are shifted in distance 
due to phase shifts at the scattering and absorbing atoms, and must be compared to standard 
compound XAFS to obtain the bond lengths. In the figures showing FT data, the shifted 
position of the peaks are given by vertical lines. In Fig. ^ and |]b, the FT of kx(k) for data 
collected on PBCO at the Pr K-edge and for Ce02 at the Ce K-edge are compared. FT's 
for several concentrations of Pr in YBCO:Pr for the Y, Pr, and Cu K-edges are presented 
in Fig. I, I, and |. 

By comparing the FT of the XAFS from the various edges, we can quickly ascertain the 
dominant substitution site for Pr. From the solid line in Fig. ^ (PBCO, Pr K-edge) and 
the solid line in Fig. ^ (YBCO, Y K-edge) we can see by the location of the first two peaks 
that the oxygen (first peak) and copper (second peak) environment around Pr in PBCO is 
very similar to that around Y in YBCO. In contrast, a simulation using FEFF5I23 of the Ba 
K-edge XAFS in YBCO puts the oxygen and copper peaks at very different positions than 
observed in the Pr K-edge data. Similarly, the Cu K-edge data for YBCO (Fig. |]) also has 
a very different environment. This indicates that most of the Pr substitutes at the Y site. 

We can also compare the Pr K-edge XAFS for the first Pr-0 peak to the Ce K-edge 
XAFS for the first Ce-0 peak in Ce02. Both sites should be 8-fold coordinated, and since 
Ce and Pr are neighbors in the periodic table, their XAFS amplitudes and phases should 
be comparable. However, the FT of kx(k) shows a 60% drop in amplitude in the Pr-0 
peak in the 100% Pr sample when compared to the Ce-0 peak (Fig. ^. In fact, the Pr-0 
peak amplitude decreases with increasing concentration (Fig. ^. The Pr-Cu peak in PBCO 
shows httle change with concentration, suggesting that the praseodymium is occupying a 
well ordered site in the unit cell while the oxygen is not. 

Changes in the oxygen environment around the yttrium in the Y K-edge XAFS, if any, 
are below the resolution of the experiment (Fig. H). Changes in the further neighbor peaks 
are only evident in the Y-Ba, Y-Y, Y-Pr regime where the double hump at 3.6 A becomes 
less lopsided toward the lower side with increasing Pr concentration. There is no visible 
changes in the Y-Cu peak. 

The Cu K-edge XAFS of PBCO show a 10% reduction in the oxygen peak amplitude 
compared to YBCO indicating that some disorder in the oxygen environment exists (Fig. ^. 
Changes in the Cu-Ba, Cu-Y, Cu-Pr multi-peak (between 2.6 A and 3.4 A) are complicated. 



with the large peak at 3.15 A breaking into two humps at 3.0 A and 3.3 A. The peak due to 
the [Cu{l)-Cu{2)]c-axis and the ab plane Cu-Cu bonds (centered at about 3.65 A) decreases 
monotonically with Pr concentration. This multi-peak includes the forward scattering due 
to the co-linear or nearly co-linear oxygens along most of these copper-copper paths. Since 
the forward scattering amplitude is strongly dependent on small deviations from co-linearity, 
the changes in this peak probably indicate small displacements of the O atoms. 

B. Background features in high-energy K-edge XAFS 

A preliminary FT of the Pr K-edge data which uses a simple spline or polynomial for fio 
indicates an anomalous hump near 1 A. By fitting the data above 6 A~^ to the Pr-0 and 
Pr-Cu peaks, and then subtracting the extrapolated fit from the original data, clear features 
in /io(E) emerge which cannot be reconciled as an XAFS signal (Fig. |^). The position of 
the feature roughly 115 eV above the absorption edge is consistent with a multielectron 
excitation, as estimated from the Z+1 model, i.e., with the N/yy transition of Nd. As 
pointed out by Holland et al.f^ a Ramsauer-Townsend-like effect can also lead to significant 
structure in the atomic background absorption. Such background features arise when the 
wavelength of the photoelectron is comparable to the dimension of an "embedded" absorbing 
atom in a solid and are essentially an XAFS of atomic origin. The effect appears to be 
larger in atoms of higher atomic number. Calculations by J. Rehr (private communication) 
show that this structure can occur at energies comparable to those of possible multielectron 
excitations. The background structure for Pr K-edge data is indeed much larger than features 
we have seen in previous work.E3 We have removed the background structure using the 
iterative procedure discussed in Ref. ^. Similar features were removed from the Ce K-edge 
data for the standard compound Ce02. Our work on these background features will be the 
subject of a future paper. 

IV. DETAILED ANALYSIS AND FIT RESULTS 

A. Fitting procedures and constraints 

All fits to the XAFS data are carried out in r-space on the FT of kx(k). We typically 
choose an appropriate range in k for the transform, and only fit it in some desired range in r. 
In the r-space fitting procedure, we vary the XAFS from standard compounds in amplitude, 
a Debye- Waller-like broadening a, position, and the K-edge threshold energy Eq to obtain 
a good fit. We have developed experimental standards for many pairs of atoms. For pairs 
which we have no experimental standards, the theoretical standards calculated by FEFF5IHZI 
were used. Fits to the Pr K-edge data used the theoretical standards, while Cu K-edge 
data were fit with experimental standards. Since we did not obtain standard compounds 
containing Pr, we evaluated the quality of the theoretical standards for Pr by comparing 
the Ce-0 standard generated by FEFF5 to experimental data obtained from Ce02, with 
good agreement for the Ce-0 peak (Fig. ^). The non-linear fitting routine minimizes a 
fit parameter, C"^ = J2 |FT(kx(k)) - fitp/ < |FT(kx(k))p >, which is roughly proportional 
to the statistical x^- Each fit can include some constraints, such as setting the number 



of neighbors in one peak equal to some fraction of the number of neighbors in another 
peak. Such constraints are necessary to restrict the number of independent parameterg22l 
and can be used to test a model, or to control correlations between fit parameters. For 
instance, correlations between amplitude and broadening factors can be controlled by, say, 
constraining the sum of the number of neighbors for two peaks to be some fixed value. The 
constraints are generally self-consistent, i.e., various atomic pair distances and number of 
neighbors are related and must have the same results for measurements at different K-edges. 

B. Determining substitution site of Pr 

A fit to FEFF5 theoretical XAFS standards that includes a mix of Pr at the Y site and 
at the Ba site indicates that, at most, the Pr exists on the Ba site 8% of the time. This 
non-zero result could only be obtained if the fit allowed the energy of the Pr K-edge (Eo) 
to be significantly different (10 eV) for Pr in the Y site and Pr in the Ba site. Such a small 
concentration of Pr at the Ba site in these fits is consistent with all Pr substituting onto the 
Y site. More information regarding this issue from the other K-edges is difficult to obtain, 
because the backscattering amplitudes of both the Pr (Z=59) and Ba (Z=56) atoms are very 
similar. 

A similar fit that allows some Pr to reside on the Cu(2) site is not as simple to interpret. 
For this discussion, we will refer to Pr at the Cu(2) site as Prc„(2), and Pr at the Y site as 
Pry. We can obtain a fit to the Pr K-edge data that includes 20% of the Pr residing on 
the Cu(2) site, with a lengthening of the Ptcu{2)-0 bonds by 0.25 A and of the Ptcu{2)-P'^y 
bonds by 0.13 A. Such a distortion in the oxygen environment is consistent with both the 
Pr K-edge and the Cu K-edge data. However, a distortion in the Pycu{2)-Pt^y bond should 
effect how the Pry sits in the unit cell, and should thus be refiected in the Pry-Cu bond. 
This effect should be large even for small (~ 20%) Ptcu{2) concentrations, because the shift 
in the Ptcu{2)-Pt^y bond length is large enough to cause nearly a complete cancellation of 
the distorted Pry-Cu peaks. The data (Fig. ^ and the fits (see Sec. |1V D|) show very little 



change in the Pr-Cu peak. Therefore, these data are not consistent with Pr at the Cu(2) site 
within the amplitude limits of XAFS spectra (~ 10%), for all Pr concentrations measured. 

C. Y K-edge fits 

Fits to the Y K-edge data were performed, without any constraints, on the 0, 30, and 
50% Pr samples. Results for the oxygen and copper peaks are in Table |I|. Further peaks, 
such as the Y-Ba, Y-Pr, and the Y-Y peak were fit but are not reported, because they are 
all at similar distances, and the Y-Ba, Y-Pr peaks are almost exactly out of phase with 
the Y-Y peak. This situation allows changes in amplitude to be mimicked by changes in 
position, and thus does not give rehable results. 

The number of oxygen neighbors shown in Table || changes very little with concentration 
(about 5% up to 50% Pr), while the copper peak shows a slight increase. Such small changes 
are not within the resolution of this experiment, and the data are therefore consistent with 
no change in the number of O atoms surrounding Y. The Y-0 distance remains constant. 



The Y-Cu distance does increase monotonically with Pr concentration, as one would expect 
from the diffraction results (Table |). 

D. Pr K-edge fits to a variety of models 

Fits to the Pr K-edge XAFS spectra were performed with and without a variety of 
constraining models and initial conditions. Each model has an associated C^ value for its 
fit, which should give the reader a feel for the range of models that may give plausible results. 



Table |ITT| gives fits to the 100% Pr sample with a variety of constraints imposed on the fits 
that simulate different models. None of the models constrain the Pr-Cu peak, and none of 
the fits to these models show any significant differences in this peak. 

Model 1 assumes a single Pr-0 bond, with no constraints, and gives a good quality of 
fit. This fit shows an oxygen deficiency around Pr, with only 4.8 of the expected 8 oxygens. 
Since this result is inconsistent with the fit results of the Y K-edge (Sec. |1V (J] ), the Cu 
K-edge (Sec. |IVE|) and the volumetric measurements of the oxygen concentration in these 
samples (Table |ip, we constrained the next fit to give 8 oxygens. Model 2, which has one Pr- 
O distance and does not allow for an oxygen deficiency, yields a poor fit with a substantially 
larger a and an increase in the goodness-of-fit parameter (C^) by more than an order of 
magnitude. It is therefore highly unlikely that simple harmonic disorder of the Pr-0 bond 
can explain the reduction in the Pr-0 peak in the FT of kx(k). 

Next we considered several models that allow the first peak to be composed of two 
different Pr-0 harmonic distributions. Since for K-edge data any shift due to differences in 
valence should be small, all of these models constrain any Eo shifts in the Pr-0 bonds to 
be the same. Model 3 allows for a double Pr-0 bond, with no further constraints. This fit 
has the best C^ of any of the models. It indicates a total number of oxygens which is still 
low (5.7). However, a fit to model 4, which constrains the total number of nearest-neighbor 
oxygens to sum to 8, gives an C^ which is only 30% higher. 

Two other models were tested to see if the fits in models 3 and 4 can be distinguished 
from similar models. All these fits have a good quality of fit parameter, but shift certain 
parameters outside a reasonable range. A fit which constrains the length of the two bonds 
to be equal and the total number of oxygens to be 8 (model 5), but allows for different cr's 
gives a very large a of 0.287 A for the short Pr-0 bond, with 3.5 O neighbors (roughly 36% 
of the O). However, a peak this broad has a tiny XAFS amplitude, and therefore this fit is 
essentially the same as the single Pr-0 peak fit in model 1. Consequently we do not consider 
this possibility further. Constraining the a's to be equal (rather than the bond lengths), 
with 8 total oxygen neighbors (model 6) gives a fit with a short Pr-0 bond length of 2.15 
A and an increase in C^ by less than a factor of two over model 4. In this fit, only 15% 
of the Pr-0 bonds are short. Such a short bond length is not unreasonable, given that the 
short Pr-0 bond in BaPrOs is 2.18 Afn although Pr is only 6-fold coordinated with respect 
to oxygen in that compound. Models 3, 4, and 6 all indicate that the radial distribution 
of Pr-0 bond lengths is not a symmetric distribution; the major weight is near 2.45 A but 
there is clearly some weight in a broad peak centered near 2.27l;o'j'2 A. 

Fits that allow for three Pr-0 distributions (as suggested by Ref. ITH) are not conclusive. 



because the number of parameters required for such a fit exceeds the maximum number of 



parameters allowed by our fitting range.E2l However, the two peak fits described above giving 
large a's for the shorter Pr-0 distribution may mimic a situation with three distances. 

The peaks that allow for two Pr-0 distributions give the most consistent results. There- 
fore, fits that allow the Pr-0 peak to be composed of two different harmonic distributions 
(models 3 through 6) suggest 15-40% of the nearest-neighbor Pr-0 bonds are either disor- 
dered and/or shifted (^ 0.18 A) compared to a more populated, ordered site. 

Table |V| summarizes the results of a two-peak fit with the constraint that there be 8 
oxygens in the first peak for all measured Pr concentrations. Debye- Waller factors for the 
disordered site increase linearly with Pr concentration, and are as much as 150% higher than 
for the "normal" site in PBCO. The fits show no appreciable change in the Cu environment 
around praseodymium with concentration, apart from the expected lengthening of the Pr-Cu 
bond due to the slight expansion of the lattice (Table |I|). The number of copper neighbors 
around the praseodymium remains the same for the 100% Pr sample compared to the other 
Pr concentration samples. Therefore, we can conclude that the effect of the small fraction of 
other phases in the PBCO sampleE3 on the XAFS is below the resolution of the experiment. 

E. Cu K-edge fits 

Cu K-edge fits to the 0% and 100% Pr samples were also carried out. Fits to partial 
concentrations of Pr are difficult to obtain because of the known distortions caused by the two 
sub-lattices. Several constraints were placed on the fits to maintain certain symmetries in the 
system. For instance, the Cu(l)-0(4) distance plus the Cu(2)-0(4) distance is constrained to 
equal the Cu(l)-Cu(2) distance obtained from diffraction.c^l The Ba position is constrained to 
be in the center of the ab-plane of the unit cell, while the Y or Pr position is constrained to be 
in the middle of all planes in the unit cell, consistent with the Cu(2)-Y and the Cu(2)-Cu(2) 
(c-axis) bond lengths. Other constraints on the number of neighbors have been imposed 
to maintain the number of atoms at certain sites while still allowing for overall shifts in 
amplitude. The proximity of the Cu(l)-0(1), Cu(2)-0(2), and Cu(2)-0(3) peaks makes an 
accurate determination of their relative amplitudes difficult. Consequently, measurements 
of changes in their relative amplitudes is beyond the resolution of our fits. We have treated 
the Cu(l)-0(1), Cu(2)-0(2), and the Cu(2)-0(3) as one peak (the Cu-Op;a„ar peak), and 
constrained the relative amplitudes of this peak to the Cu(l)-0(4) and the Cu(2)-0(4) peaks. 

The fits to the Cu K-edge data are summarized in Table 0. Fit parameters for YBCO 
are typical of fits we have obtained in previous work. There are several differences between 
the YBCO and the PBCO results. Bond length changes are, of course, expected, and will be 
dealt with in some detail in the next section. Number of neighbor values are 15% different 
for the combined Ba and Pr peak compared to YBCO, which is usually considered to be 
outside the expected error (10%) for XAFS amplitudes. However, the amplitude of this 
combined peak is very sensitive to the individual positions of the Ba and the Pr, as well 
as to the position of the Cu(2)-Cu(2)c_axis peak, and are therefore not very reliable. The 
broadening factors are slightly larger for the planar oxygens and the Cu(l)-0(4) bond in 
PBCO, as are the cr's for the Ba and Pr peaks. The Cu(2)-0(4) bond in PBCO is much 
narrower than in the YBCO. Again, the a's for the Ba and Pr peak suffer from the same 
problems as their amplitudes, but the change in a for the oxygen peaks is probably real. 
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V. DISCUSSION 
A. Consistency bet^veen various K-edge data 

Since each absorption edge gives us information about the oxygen environment within 
YBCO:Pr, we can constrain our results to give a consistent picture. Y K-edge data show 
httle if any change in the structure of YBCO:Pr up to 50% Pr, except (possibly) in the 
expected lengthening of the Y-Cu bond. Even this expansion of the Y-Cu distance is at the 
edge of the accuracy of this experiment. These data therefore show that up to 50% Pr, the 
oxygen and copper environment around yttrium in YBCO: Pr is essentially the same as in 
pure YBCO. More specifically, there is little evidence of a loss of oxygens or of disorder in 
the oxygen or copper environment. 

The Pr K-edge data and fit results indicate that either there are missing oxygens in the 
Cu-0 planes or that there is a large amount of disorder and/or distortion in these planes. 
The former is inconsistent both with the measured oxygen content and with the number of 
neighbors for the Cu-0 and the Y-0 peaks in the XAFS data, and is therefore ruled out. 
The data is best fit by two Pr-0 distributions, separated in pure PBCO by about 0.18 A, 
with one distribution less populated and more disordered than the other. 

The Cu K-edge data is more complicated because of the two copper sites. Nevertheless, 
the oxygen environment as seen from the coppers should be consistent with (and therefore 
can help constrain) the Pr K-edge results. The fit results to the oxygen peaks show essentially 
no change in the number of nearest-neighbors, but indicate a broadening in the Cu-Opianar 
peak consistent with the 10% decrease in the peak height. If the Pr-0 distortion is completely 
radial, then a shortening of 0.18 A would correspond to only a 0.04 A lengthening of ~30% 
of the Cu(2)-0(2,3) bonds, which is not resolvable in the fits. However, such a distortion will 
contribute to a broadening of the Cu- Opianar peak, and therefore the measured broadening 
is consistent with the proposed radial distortion in the Pr-0 bond. The small change in 
the number of Y-0 neighbors (which can also be modeled as a slight increase in a) is also 
consistent with a radial distortion in the Pr-0 bonds. 

B. On the question of valence 

The fit results to the first oxygen peak from all three edges are consistent with a distortion 
of the Pr-0 bond in the radial direction of approximately 0.18 A. This distortion may be 
due simply to a hybridization of the O 2p and Pr 4/ electrons. However, the magnitude 
of the shift is very consistent with some Pr^"''; the best fit gives a Pr-O^hort bond length 
(2.27 A) which is close to the nearest-neighbor bond length in 8-fold coordinated Pr02. By 
constraining the broadening factors, we can obtain a much shorter bond length (2.15 A) 
which is more consistent with 6-fold coordinated BaPrOs. The comparison to the 8- fold 
coordinated Pr02 makes the most sense, since Pr in PBCO is 8-fold coordinated, although 
it is possible that not all the oxygens participate in Pr-0 bonds. 

If we assign all of the short Pr-0 bond lengths to the formal Pr^+ state, our fit results 
indicate that the Pr is in that state 15-40% of the time, corresponding to a formal valence of 
-|-3.33io;?g. Such an assignment is reasonable for the shortest measured Pr-0 bond length. 
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but if the bond length is closer to ~2.3 A, then some could still include Pr^"*", as in bixbyite 
Pr203. Therefore, the upper limit on the valence range may be too high. Likewise, the lower 
limit was obtained by forcing the broadening factor of the short bond to be much lower than 
in the best fit, which may not be reasonable. 

The possible disorder in the shorter Pr-0 bond could cause localization of charge carriers, 
and could be a byproduct of a 4/-2p hybridization. On the other hand, since the disorder 
is greatest for the PBCO sample, the disorder could be due to the less common Pr^^-Pr^"*" 
combination, which will probably have different Pr-0 bond lengths than either the Pr^+-Pr3+ 
or the Pr3+-Pr^+ combinations. 

Our measured dominant Pr-0 distance of 2.45±0.01 A is in excellent agreement with 
diffraction resultgi2rH3 that give a mean Pr-0 distance of 2.4539±0.0038 A. We see no 
evidence of a longer distance (~ 2.5 A) as suggested by the Pr-0 bond length trend with 
ionic radius.E3. Fits that include such a peak invariably shift the peak back to the values 
we report. Therefore, we must conclude that the discrepancy of this bond length with ionic 
radius is not directly due to a mixed valency. Indeed, given the existence of a short Pr-0 
bond, the long Pr-0 bond seems to be entirely Pr^+(4/)-0. 

The range for the valence is in approximate agreement with the electronic studies men- 
tioned in Sec. | that give the valence of Pr to be close to 3-I-. However our results clearly 
suggest a higher formal valence than 3.0, as was found in Pr L-edge studies of Horn et al. 
(3.1+ for PBCO)§ and Lytle et al. (3.45+ for 20% Pr and 3.25+ for 60% Pr).0 In addi- 
tion, our structural results show that the number of short Pr-0 bonds, that is, formal Pr^+, 
increases with increasing Pr concentration, and that any Pr^+ that may be present resides 
on the Y site. Both of these results are in contrast to the Lytle et al. results. 

Our result is in striking agreement with the 4+ / hybridization model put forth by 
Fehrenbacher and Ricejla which predicts that the Pr^+ state will exist 60-70% of the time. 
A ligand hole localized on the Pr-0 (2) /O (3) bond would remove a hole from the conduction 
band and yield a formal Pr^"*" site. The resulting Pr^^-O bond length would likely be 
significantly shorter than the Pr^^-O bond, based on the observed bond lengths in other 
materials. Although we cannot give a precise estimate of the number of short bonds present, 
our results are consistent with their prediction. 

C. On the question of the relative position of 0(4) 

The Cu K-edge data for PBCO are consistent with the diffraction results of Lowe-Ma 
and VanderahBil which show the axial oxygen (0(4)) to be pushed toward the planar Cu site 
when compared to the YBCO structure (Table P^ ). The comparison to YBCO is properly 
made by taking the differences in the lattice parameters a, b, and c into account. Therefore, 
we also show the bond-lengths for "pseudo-YBCO" which is calculated using the lattice 
parameters for PBCO from Lowe-Ma and VanderahEl and the relative atomic positions of 
YBCO from Beno et al^ The difference in Cu(2)-0(4) distance between the diffraction and 
the XAFS results can be explained by comparing the Cu(2)-0(4) distance as a function 
of oxygen concentration for R123-type materials. In all these materials that have been 
measured, the Cu(2)-0(4) distance decreases with increasing oxygen concentration. Our 
result is consistent with this trend (see Ref. |2T]). An extrapolation of the measured Cu(2)- 
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0(4) distance as a function of oxygen concentration in PBCOjHll to an O content of 6.98, 
yields a bond length of 2.23 A, in excellent agreement with our result of 2.22 A. 

D. On the question of clustering in RBCO:Pr (R=rare earth and Y) 

An interesting property of Ri_a;Pra;Ba2Cu307 is that the larger the rare-earth radius, 
the lower the Pr concentration required to completely suppress Tc. It has been suggestedc3 
that this is evidence for clustering of Pr in YBCO:Pr (where the ionic-size difference is 
comparatively large), thereby maximizing the superconducting Y-Y pairs of cells within a 
crystal. We had expected that both the Pr and Y K-edge data would be very sensitive to 
this clustering, by giving the number of neighboring Y and Pr atoms. The data in Fig. ^ 
and ^ show a significant change in the structure near 3.5 A where the Y or Pr neighbor peak 
occurs, compared to the pure materials YBCO and PBCO. This means that in the 30 and 
50% samples there are both Y and Pr neighbors at each Y/Pr site. If large scale clustering 
were present, these significant changes in the XAFS would not occur. Unfortunately, fits 
that try to precisely determine the ratio of Pr to Y neighbors are plagued by the interference 
effects from the Y-Ba and Y-Y peaks, as described in Sec. |IVC|, which reduces our sensitivity. 



Consequently, we cannot differentiate between a random distribution and small clusters. 



VI. CONCLUSION 

In summary, we have collected XAFS data on samples with several different concentra- 
tions of Pr in Yi_a;Pra;Ba2Cu307 and on Ce02 at the Y, Pr, Cu, and Ce K-edges. The main 
result of this investigation is that the XAFS data and analysis show that disorder in the 
oxygen environment around the praseodymium clearly exists. This disorder is most obvious 
when comparing the first oxygen peak in Ce02 to the first oxygen peak in PrBa2Cu307, 
which shows a 60% reduction in the Pr-0 peak. The amplitude of this peak clearly shrinks 
with increasing Pr concentration. 

This PBCO data is best fit by two harmonic radial distributions for the Pr-0 bond, 
which are centered at 2.27tai2 A and 2.45±0.01 A in PBCO. The shorter bond may be 
broader than the longer one by as much as a factor of 3. The magnitude of this distortion 
agrees well with the change in bond length in other Pr oxides between Pr'^^-O (2.18 A - 2.32 
A) and Pr^+-0 (2.33 A - 2.66 A) and is therefore the first structural evidence of a split in the 
Pr-0 bond lengths, and strong evidence for a mixed valent state of Pr in Yi_a;Pra;Ba2Cu307. 

We have also verified some results of previous measurements. By comparing the Pr K- 
edge data to FEFF5 simulations of PrBa2Cu307 with Pr at the Y site and at the Ba or 
Cu(2) site, we see very little evidence for any of the Pr existing on the Ba or the Cu(2) site 
(upper limits of 8 and 10%, respectively). In addition, Cu K-edge data indicate the 0(4) site 
has moved closer to the Cu(2) site (in agreement with Ref. |21|) and therefore may impede 
charge transfer between the planes and the chains. 
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FIGURES 

FIG. 1. Tc vs. Pr concentration. The inset shows the structure of YBa2Cu307 with the 
notation used in this paper. 

FIG. 2. kx(k) for the 50% Pr sample from the Y, Pr and Cu K-edges (at 17080 eV, 41991 eV, 
and 8979 eV, respectively). 

FIG. 3. Fourier Transforms of kx(k) for (a) PrBa2Cu307 and (b) Ce02, together with fits 
to the first two major peaks (dotted). The outer envelope is the magnitude of the FT, while 
the modulating curve is the real part. All the fits discussed in the text are of similar quality 
to the PBCO fit shown here. The first peak in both cases is an eight-fold coordinated oxygen 
peak. Single pair peak positions are approximately given by vertical lines that meet the x-axis. 
These peak positions are shifted from the actual positions by phase shifts at the absorbing and 
backscattering atoms. Both transforms are from 3.5-17.0 A~^, with a 0.3 A~^ gaussian window. 

FIG. 4. Fourier Transform of kx(k) for Y K-edge data for 0% (sohd), 30% (dotted), and 50% 
(dash) Pr concentrations. The only dramatic changes occur in the Y-Ba, Y-Y, Y-Pr region. These 
transforms are taken from 3.5-15.5 A~^, with a 0.1 A~^ gaussian window. 

FIG. 5. Fourier Transforms of kx(k) for Pr K-edge data for 30% (solid), 50% (dotted), and 
100% (dash) Pr concentrations. Notice the decrease of the first peak with Pr concentration. This 
decrease is most likely due to a combination of the relative amplitudes of two Pr-0 distances and 
additional disorder in the shorter bond. The second peak is the Pr-Cu peak. It shows relatively 
little change with Pr concentration, indicating that the Pr is well ordered with respect to the Cu, 
and thus to the unit cell. The peaks in the 3-4 A range are due to a mix of Pr-Y, Pr-Pr and 
Pr-Ba. This region is difficult to fit accurately because the Y and Pr backscattering amplitudes 
are similar, but vr out of phase. The FT ranges are from 3.5-17.0 A~^, with a 0.3 A~^ gaussian 
window. 

FIG. 6. Fourier Transforms of kx(k) for Cu K-edge data for 0% (sohd), 30% (dotted), 50% 
(short dashed), and 100% (long dashed) Pr concentrations. The first Cu-0 peak decreases with 
increasing Pr concentration by about 10% in PBCO. The Cu-Y, Cu-Pr, Cu-Ba region also changes 
dramatically, due to the changing phase shifts between these peaks; see text for further discussion. 
These transforms are taken from 3.0-15.8 A~^, with a 0.1 A~^ gaussian window. 

FIG. 7. "Free-atom absorption" fio (solid) and total absorption fi (dotted) coefficients times 
sample thickness, t, as a function of energy above the K-edge. /Xot was determined by an iterative 
procedure.E2l 
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TABLES 



TABLE L Sample information: a, b, and c lattice parameters for the Yi_3;Pra;Ba2Cu30y sam- 



ples. 



y 



0.0 
0.3 
0.5 
LO 



3.821 
3.826 
3.846 
3.902 



3.882 
3.894 
3.907 
3.916 



11.621 
11.667 
11.691 
11.715 



6.96 
6.94 
6.96 
6.98 



TABLE IL Y K-edge fit results. These fits include further peaks than the Y-Cu peak, but these 



peaks are plagued by an interference effect when Pr is present (see Sec. |VD|). Y-Cu amplitudes 
are consistently larger than those calculated by FEFF5 by nearly 25%. Since we do not have data 
for a similar standard pair, we have normalized the Y-Cu number of neighbors to the YBCO result. 
The number of Y-0 nbrs is not normalized. The fits are from 1.0-3.1 A in R, and from 3.5-17 A~^ 
in the wave vector k. Each bond is expected to be comprised of 8 neighbors. The errors are ±0.01 
A in R, and ±10% in both a and number of neighbors for all XAFS fits reported in this work, 
unless otherwise noted. 







Y-0 






Y-Cu 




X 


R(A) 


nbrs 


^(A) 


R(A) 


nbrs 


a{k) 


0.0 


2.41 


7.3 


0.058 


3.21 


8.0 


0.058 


0.3 


2.41 


7.1 


0.056 


3.22 


8.6 


0.060 


0.5 


2.41 


6.9 


0.057 


3.23 


8.6 


0.060 



TABLE III. Fit results on PrBa2Cu307 to a variety of models. Fitting procedures are detailed 
in the text. C^ is the fitting parameter, and is approximately proportional to the statistical x^i 
as described in the text. The r-space fits are from 1.0 - 3.0 A in R, and the FT from 3.5-17.0 
A~^ in wave vector k. The number of Pr-0 neighbors is normalized to the Ce02 result. No 
normalization has been applied to the Pr-Cu peak. Model 1 only allows for one Pr-0 distance, 
with no constraints on the fit parameters. Model 2 also only allows for a single Pr-0 bond, but 
constrains the amplitude to be 8, as is expected from the known crystal structure. Model 3 allows 
for two Pr-0 bonds (Pic-Oshort and Pr-O^o^^, as do all the subsequent models) with no further 
constraints. Model 4 constrains the total number of Pr-0 bonds to 8. Model 5 constrains the 
number of Pr-0 bonds to 8 and holds their bond lengths equal. Model 6 constrains the total 
number of Pr-0 bonds to 8 and holds their broadening factors (a) equal. 





P] 


"-Oshort 






Pr-Ozong 






Pr-Cu 






model 


R(A) 


nbrs 


a{k) 


R(A) 


nbrs 


ct(A) 


R(A) 


nbrs 


a{k) 


C2 


1 








2.46 


4.8 


0.066 


3.27 


7.8 


0.049 


0.148 


2 








2.47 


8.0 


0.116 


3.27 


7.6 


0.047 


1.513 


3 


2.30 


1.0 


0.135 


2.46 


4.7 


0.063 


3.27 


7.9 


0.048 


0.132 


4 


2.27 


2.9 


0.168 


2.45 


5.1 


0.067 


3.27 


7.7 


0.048 


0.174 


5 


2.46 


3.5 


0.287 


2.46 


4.5 


0.064 


3.27 


7.8 


0.048 


0.160 


6 


2.15 


1.2 


0.084 


2.44 


6.8 


0.084 


3.27 


7.7 


0.048 


0.300 



TABLE IV. Pr K-edge fit results to Yi_:rPrxBa2Cu307 as in model 4 in Table |T|. 



X 


R(A) 


Pr-Oshort 

nbrs 


a(A) 


R(A) 


Pr-O^ong 

nbrs 


a(A) 


R(A) 


Pr-Cu 
nbrs 


a(A) 


0.3 
0.5 
1.0 


2.25 
2.29 
2.27 


1.7 
2.7 
2.9 


0.128 
0.115 
0.168 


2.43 
2.43 
2.45 


6.3 
5.3 

5.1 


0.056 
0.053 
0.067 


3.25 
3.25 
3.27 


7.5 
7.4 

7.7 


0.046 
0.047 
0.048 
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TABLE V. Cu K-edge fit results to YBa2Cu307 (YBCO) and PrBa2Cu307 (PBCO). The 
"pseudo-YBCO" is calculated from the lattice parameters of PBCOEll and the relative atomic 
positions in YBCO.I^ The last column shows the differences between the XAFS results and the 
diffraction results for PBCO. The samples used in the diffraction measurements for PBCO had a 
mean oxygen content of ~ 6.77. The PBCO XAFS results agree with the diffraction results, within 
calculated errors, except for the Cu(l)-0(4) and the Cu(2)-0(4) bond lengths, as discussed in Sec. 
^. The calculated errors are the same as discussed in Table || except for the Cu-Ba,Pr neighbors 
in PBCO (± 1.5), and the positions of the Cu-Ba, Pr and Y peaks for both PBCO (± 0.03 A) and 
YBCO (± 0.02 A). 

expected PBCO XAFS YBCO XAFS pseudo-YBCO PBCO diffrac.0 XAFS-diffrac. 

bond nbrs R(A) nbrs a{A) R(A) nbrs o-(A) R(A) R(A) R(A) 

1.849 +0.031 

1.965 +0.015 

2.254 -0.032 

3.388 -0.018 

3.482 -0.012 

3.265 -0.025 



Cu(l)-0(4) 


2 


1.88 


2.4 


0.082 


1.85 


2.3 


0.069 


1.856 


Cu-Opianar 


10 


1.98 


11.8 


0.082 


1.94 


11.7 


0.069 


1.965 


Cu(2)-0(4) 


2 


2.22 


2.4 


0.056 


2.27 


2.3 


0.110 


2.311 


Cu(2)-Ba 


8 


3.36 


9.5 


0.056 


3.37 


7.6 


0.035 


3.408 


Cu(l)-Ba 


8 


3.47 


9.5 


0.056 


3.49 


7.6 


0.035 


3.499 


Cu(2)-Y,Pr 


8 


3.23 


9.5 


0.051 


3.21 


7.6 


0.051 


3.232 
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